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CHAPTER 1: INTRODUCTION 


There has been a rapid growth in the use of polymers and polymeric-based materials in 
recent years, and there is every indication that this trend will continue. Of particular interest has 
been the use of polymers as the matrix component in structural composite materials. Polymers 
may be roughly categorized as being either thermosets or themoplastics. Both thermosets and 
thermoplastics are comprised of long molecular chains. However, the molecular chains of 
thermosets are highly cross-linked, forming an extensive three-dimensional molecular structure. 
Conversely, thermoplastics do not form cross-links. One of the ramifications of this difference in 
molecular structure is that once a thermoset polymer has been polymerized (i.e., once the cross- 
links have been formed) the thermoset cannot be melted. On the other hand thermoplastics can be 
readily melted and remolded by the application of heat and pressure. Most thermoplastics have 
glass-transition and melting temperatures much lower than the glass-transition temperature of 
thermosets, and thus thennoset resins have been used most often in structural composite materials. 
However, recent advances in polymer materials technology have resulted in a new generation of 
high-temperature thermoplastic polymers which are better suited for engineering applications. 
Examples of these new-generation thermoplastics include poly-ether-ether-ketone (“PEEK”, with a 
melting temperature of roughly 340°C) or a variety of polyimide resins (with a melting temperature 
ranging from 230-350°C). Because of these thermal properties as well as the inherent toughness 
and relative ease with which thermoplastics can be fabricated, recent research and developmental 
efforts in structural composite materials have turned to thermoplastic-based matrices. 

The expanding role of thermoplastic polymers in load-bearing structural applications has 
prompted a need to develop a design methodology to predict the mechanical response of these 
materials under complex loading conditions. In particular, prediction of yield and post-yield 
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behavior is of current technical interest. At present most of the design tools and concepts used in 
the analysis of polymeric- based materials are based on the same principles used for metals and 
metallic alloys. However, the atomic/molecular structure of polymers is entirely different from that 
of metals. Metals are chjiracterized by a highly crystalline atomic structure whereas polymers are 
characterized by long molecular chains that may be amorphous or semi-crystalline at the molecular 
level. Since the macroscopic yield and post-yield behavior of any material is fundamentally 
governed by its atomic/molecular microstructure, it is unlikely that a design methodology used for 
metals will apply equally well to polymers. Therefore, it may be necessary to develop new design 
methodologies for polymers, or to modify existing design tools that had been developed for use 
with metals. 

Purpose 

The purpose of this research is to study the yield and post-yield behavior of a thermoplastic 
polymer subjected to biaxial stress states. In particular, experimental results will be compared with 
theoretical predictions based on classical plasticity theories previously applied to metals. Of 
particular interest is whether the yield or post-yield behavior of thermoplastic polymers is 
influenced by the hydrostatic stress component 

The material chosen for study was high-density polyethylene. High-density polyethylene 
has a melting temperature of roughly 130°C, and is not considered to be a potential matrix material 
for use in structural composites. Polyethylene was selected because (1) it is a semi-crystalline 
thermoplastic, and hence its post-yield behavior may be representative of the general class of 
thermoplastic polymers, (2) it is readily available in the form of thin- walled tubes, and (3) it is very 
inexpensive as compared to new-generation thermoplastics such as PEEK or polyimides, and 
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hence a large number of tests could be conducted at modest cost 

The study involved plane stress problems where thin-walled tubes of high-density 
polyethylene were tested under combined tensile axial loads and internal pressures at room 
temperatures. A rather specialized testing apparatus was developed during the study in order to 
perform these tests. This testing system is now fully calibrated and operational, and will be used 
to study the inelastic behavior of other new-generation thermoplastic resins of greater structural 
interest during the coming year. 



CHAPTER 2: BACKGROUND INFORMATION 


THE PREDICTION OF YIELDING 
Definition of the Yield Stress 

Conceptually, the yield point of a material defines the transition from purely elastic to 
elastic/plastic behavior. F or some materials, the yield point is indicated by a sharp drop in stress in 
a uniaxial stress-strain curve, followed by slowly increasing stress with increasing strain. In these 
cases the near-discontinuity in stress makes identification of the yield stress very easy. For other 
materials no such stress drop occurs, but instead the slope of the stress-strain curve gradually 
decreases from an initial "elastic" value, which defines Young's modulus, to a lower value. In 
these cases the identification of yielding becomes more problematic. A common method of 
defining yield in this case is to use a permanent strain offset. For metals, this strain offset is 
usually defined as 0.2% strain. 

Defining yield in polymers involves complexities not normally encountered with metals. 
Polymers are viscoelastic and therefore the measured yield stress is sensitive to the stress- or 
strain-rate imposed. Researchers have attempted to define the yield stress of polymers in a variety 
of ways. One method [Bowden and Jukes, 1972; Freire and Riley, 1980; Mears et al, 1969; 
Whitfield and Smith, 1972] is identification of the point of load drop, as described above. Another 
method [Carapellucci and Yee, 1986] is to perform a constant stress creep test, and to define yield 
based upon the intersection of tangents drawn to the initial and final portions of the resulting strain 
versus time curve. Probably the most common method is to apply a monotonically increasing axial 
load at a constant rate, and to define yielding on the basis of a strain offset. However, no standard 
loading rate or strain offset has been established. Offset strains ranging from 0.3% strain [Caddell 
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and Woodliff, 1977; Raghava et al, 1973] to as high as 2% strain [ Ely, 1967; Pae, 1977] have 
been used. 

Obviously, the different methods to determine yield result in distinctly different values of 
measured "yield strength" for a given polymer. Figure 2.1 compares various methods used to 
determine the yield strength of a polymer, and indicates that distincdy different values of the "yield 
point" can be deduced from the same experimental data. The inconsistency in the definition of 
"yielding" makes comparison of results obtained during different studies difficult. A standard 
definition, such as exists for metals, needs to be adopted for use with polymers. 

The method used to define yielding during the present study will be discussed in detail in a 
following section. At this point it is appropriate to note that yielding was defined on the basis of a 
0.3% offset in the octahedral shear stress vs octahedral shear strain curve of polyethylene. A value 
of 0.3% was selected because it is a widely used value of strain offset, at least within the polymers 
technical community. Yielding was defined in terms of octahedral shear stress and octahedral 
shear strain because this approach automatically accounts for biaxial loading effects. Also, all tests 
were conducted under a constant rate in order to minimize rate effects . Specifically, the yield and 
post-yield response of the polyethylene specimens was measured under a condition of constant 
octahedral shear stress loading. 

Yield Criteria 

The yield criteria which have been developed for metals typically make the following 
assumptions: (1) the initial compressive and tensile yield strengths are equal, (2) the hydrostatic 
stress component does not contribute towards yielding (i.e., the deviatoric stress component 



xial Stress (psi) 



Figure 2.1 - Comparison of Various Methods of Determining the Yield Point 

(Data Taken from Load Path 1) 
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governs yielding), and ( 3) the material is isotropic and homogeneous. If the third assumption is 
relaxed the resulting yield criterion is suitable for use with anisotropic materials. The most 
commonly used isotropic and anisotropic yield criterion will be briefly reviewed in seperate 
sections below. Both "types" of yield criterion involve the hydrostatic and octahedral shear stress 
levels. Hydrostatic stress (or "mean" stress) may be written in terms of the three principal stresses 
as: 


(ai + 02 + 03) 
^hyd — 'i 


( 2 . 1 ) 


The deviatoric stress level associated with some given stress state is characterized by the octahedral 
shear stress. The octahedral shear stress may be written in terms of the principal stresses as: 

Xoct = j V (ai- 02 ) 2 +(a 2 - 03 ) 2 -Kcr 3 -ai ) 2 (2.2) 

Isotropic Yield Criteria: The isotropic yield criterion most widely used for ductile metals and 
polymers is probably the Von Mises criterion. According to the Von Mises criterion, yielding 
occurs when the octahedral shear stress in the material reaches a critical value: 

(Oi - 02 ) 2 + (02 - 03)2 + (03 - O 1) 2 = 2T 2 (2.3) 

where 01 > 02 > 03 and T is the uniaxial tensile yield strength of the material. 

The Von Mises yield criterion is based on the assumptions that the hydrostatic stress 
component does not contribute to yielding and that the compressive and tensile yield strengths are 
equal. Experimental evidence indicates that many polymers violate both of these assumptions. 
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That is, polymers do in general have distinctly different compressive and tensile yield strengths and 
furthermore their yield behavior is significantly influenced by the hydrostatic stress level. For 
example, Christiansen, Itaer, and Radcliffe [1971] found that the yield strength of polycarbonate, 
polyethylene terephthalate, polychloro-trifluoroethylene, and polytetrafluoroethylene all increased 
with increasing hydrostatic pressure. Mears, Pae and Sauer [1969] have reported similar pressure- 
dependent yield behaviors for polyethylene and polypropylene. The same conclusions have been 
reached for a variety of other polymers during studies by Ainbinder et al [1965], Pae and Mears 
[1968], Sardar et al ([1968], and Spitzig [1979]. This experimental evidence indicates that yield 
criteria which assume that the hydrostatic stress component does not influence the yield behavior, 
and/or which assume that the initial compressive and tensile yield strengths are equal are not 
applicable for use with polymers. 

A modified version of the Von Mises criterion, in which these two assumptions are 
relaxed, was first presented by Stassi-DAlia [1967] and subsequently presented in a different form 
by Raghava, Caddell, and Yeh [1973]. Raghava et al modified the Von Mises criterion in the 
following manner: 

(°i - o 2 ) 2 + (02 - 03)2 + (03 - Oi ) 2 + 2(C - T)(oi + o 2 + 03 ) = 2CT (2.4) 

where C and T are the uniaxial compressive and tensile yield strengths respectively. This criterion, 
known as the Pressure-modified Von Mises criterion, is based on the assumption that yielding is a 
function of both the hydrostatic and deviatoric stress states. The first three terms are equivalent to 
the octahedral shear stress term of the original Von Mises criterion, while the last term is associated 
with the hydrostatic stress component . 
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Note that effects of hydrostatic stress and differences in compressive/tensile strength are 
not independent in the Pressure-modified Von Mises criterion. That is, if the compressive and 
tensile yield strengths of a material are equal then the hydrostatic stress term disappears and the 
Pressure-modified Von Mises criterion reduces to the original Von Mises criterion. It is interesting 
to note that a similar interdependence between tensile/compressive yield strengths and hydrostatic 
effects has been demonstrated for some metals. Spitzig, Sober, and Richmond [1976] have shown 
that some sintered powder materials (e.g., sintered iron) exhibit significant differences in 
compressive and tensile yield strength, and also that the yield behavior of these metals is influenced 
by the hydrostatic stress component. Similar conclusions were reached by Betten, Frosch, and 
Borrmann [1982] for the case of quenched and tempered steels. 

The yield behavior of polymers has been compared to both the Von Mises and Pressure- 
modified Von Mises Criterion. Polymers which have been studied include acrylics [Ely ; 1967], 
polyester mixtures [Freire and Riley ;1980], and polycarbonate, polyvinylchloride, and 
polyethylene [Raghava and Caddell, 1973; Raghava et al, 1973]. In all cases the Pressure- 
modified Von Mises criterion was a better predictor of yield than the original Von Mises criterion, 
indicating that the yield response of these polymers is sensitive to the hydrostatic stress level. 

Anisotropic Yield Criteria: Anisotropic materials possess different material properties in 
different directions. Until recently, most research on the yield behavior of materials has focused 
on isotropic models. However, many materials are anisotropic. In addition, when shaped by 
manufacturing processes such as rolling or extrusion even materials which are initially isotropic 
may become oriented at the atomic/molecular level and therefore become anisotropic at the 
macroscale. This is particularly true of polymers since they are composed of long molecular chains 
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which can easily become oriented in one direction, resulting in a polymer with distinctly different 
stiffnesses and strengths in different directions. 

The criterion used most often to model the behavior of anisotropic materials is the Tsai-Hill 
yield criterion . Hill [ 19 * 50 ] developed this criterion by selecting a form which would reduce to the 
Von Mises criterion if the yield strengths were the same in all directions. Tsai has subsequently 
applied this theory to composite materials, and the resulting theory is now known as the Tsai-Hill 
yield criterion. It can be written in the following manner in terms of the principal stresses Gj, 02, 
and 03: 


H(gi - G2) 2 + F(G2 - G3) 2 + G(G3 - Gi) 2 = 1 


with: 


H + G = 


1 

T1 2 


H + F = 


T 2 2 


G + F = -^ 
T3 2 


( 2 . 5 ) 


Since the material is anisotropic, the coefficients in the equation are functions of the tensile yield 
strengths in the 1 -, 2 - and 3 - directions, Ti, T2 and T3 respectively. If the material is transversely 
isotropic, then T2 and T3. are equal and consequently G is equal to H. The Tsai-Hill yield criterion 
is related to the maximum shear which exists in each principal plane; it includes no linear terms 
which would represent the hydrostatic stress level, and also ignores any differences in initial 
compressive and tensile yield strengths. Stassi-D'Alia [ 1969 ] modified the Tsai-Hill criterion to 
include these factors, and Caddell, Raghava, and Atkins [ 1973 ] proposed a similar yield criterion 
at a later date. The Pressure-modified Tsai-Hill criterion, as it is sometimes known, is defined as 
follows: 


H(Oi - O2) 2 + F(G2 - O3) 2 + G(G3 - Gi ) 2 + KiGj + K2G2 + K3G3 = 1 


(2.6) 
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with H + G = qtj 

H + F = 7 i- 
C2T2 

g + f -c 3 t 3 

K, C 1' T 1 

Kl “ C1T1 

Kt - C 2 " T 2 

K2_ c 2 t 2 

K - c l:Jl 
K3_ C 3 T 3 


The compressive and tensile yield strengths in the 1-, 2-, and 3- directions are denoted by Cl, Tj, 
C2> T 2 , C 3 , and T 3 . If the material is transversely isotropic then C 2 and T 2 are equal to C 3 and 
T 3 , and therefore, G = H. If the compressive and tensile strengths are equal, Eq (2.6) reduces to 
the original Tsai-Hill criterion, and if material properties are the same in all directions, it further 
reduces to the original Von Mises criterion. 

The few experimental studies which have been undertaken to investigate the yield behavior 
of anisotropic polymers seem to validate the Pressure-modified Tsai-Hill yield criterion. Studies 
by Rider and Hargreaves [1969] and Shinozaki and Groves [1973] followed methods suggested 
by Hill in which the change in the tensile and compressive yield strengths of oriented polymer 
sheets was measured as a function of the angle of orientation 0. Caddell and Woodliff [1977] 
conducted a series of tests using highly anisotropic specimens of polycarbonate, polyethylene, and 
polypropylene. The test specimens were oriented by applying a tensile load to a specimen until a 
stable neck had formed The final test specimens were then machined out of these necked regions. 
A good correlation between experiment and theory was reported. Carapellucci and Yee [1986] also 
conducted tests on anisotropic polycarbonate specimens and reported good agreement. Hence, it 
appears that the yield behavior of anisotropic polymers is significantly affected by the hydrostatic 
stress level. 



THE PREDICTION OF POST- YIELD BEHAVIOR 


Post- yield behavior typically involves large deformations, such that the familiar definitions 
of conventional "engineering" stress and strain no longer apply. Conventional engineering stress a 
and engineering strain e may be converted to true stresses s and true strains e as follows: 


s = a(e+l) (27) 

e = ln(e+l) (2.8) 

A true stress-true strain curve is often called a plastic "flow" curve because it gives the 
stress which causes the material to flow plastically for any given strain. Post-yield theories are 
typically developed in terms of the true stress and true strain, and in this report true stress and true 
strain will be used to describe the post-yield behavior of polyethylene. 


The post-yield theories considered during this study will be briefly reviewed below. Both 
isotropic and anisotropic materials will be discussed. It is appropriate to note that, in contrast to 
the yield criterion reviewed above, all existing post-yield criteria assume material incompressibility 
during plastic flow. The possibility that the flow behavior of polymers is influenced by the 
hydrostatic stress level has apparently not yet been addressed within the technical community. 
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Post Yield Behavior of Isotropic Materials 

The assumption that plastic deformations are independent of the hydrostatic stress level 
leads to the conclusion that the volume of a solid remains constant during plastic flow. Therefore, 
during plastic flow the true strains are related by: 

ei+e2+e3=0 (2.9) 

Equation (2.9) represents the "constancy of volume" assumption. Note that if v = 1/2, then a 
similar relationship holds for engineering strains: 

e 1 + £2 + £3 = ( l-2v)(a! + 02 + a 3 )/E = 0 (when v = 1/2) (2.1 0) 

For many metals undergoing plastic deformation induced by a uniaxial stress, the flow curve can 
be expressed by the simple power law relation, also known as the Ludwik expression: 

s = Syield + Me n (2.11) 

where s yield is the true stress at yield, M is the strength coefficient and n is the strain hardening 
exponent. This equation is only valid for isotropic materials and is defined from the onset of 
plastic flow to the maximum load at which the specimen begins to neck. Equation (2.1 1) provides 
a simple mathematical expression to describe the post-yield flow curves of a material if constants M 
and n are known. 
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Equation (2.11) is appropriate for use under uniaxial loading conditions. In a complex 
state of stress and strain it i s more useful to relate invariant functions of true stress and true strain. 
The invariant functions most frequently used to describe plastic deformation of metals are the 
effective true stress and effective true strain [Dieter, 1986]. The effective true stress is defined as a 
function of the principal true stresses: 

l 

s = -p [(si-S2) 2 + (S2-S3) 2 + (S3-S1) 2 ] 2 (2. 12) 

4l 

The effective true strain is defined as a function of the principal true strains: 

1 

e = -^ [(ei-e2) 2 + (e2-C3) 2 + (e3-ei) 2 ] 2 1 3) 

The power law relation caa now be expressed in the form of effective true stress and effective true 
strains: 


s = Syieid + Me n 

An assumption implied by Equation (2.14) is that the flow curve is independent of the hydrostatic 
or mean stress component. That is, the post-yield behavior is assumed to be governed exclusively 
by the deviatoric stress tensor. 

The power law expressions (i.e., either Eqs. 2.11 or 2.14) provide a simple mathematical 
expression to describe the flow curves of a material, and the power law formulation has been 
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widely applied because of this simplicity. However, these relations are entirely empirical and 
deviations from the power law are frequently observed. 

More rigorous expressions to describe the post-yield behavior of a material have been 
developed. Two such systems of equations are the Levy-Mises and Prandtl-Reuss equations. In 
either case two fundamental assumptions are made: it is assumed that the material is elastic - 
perfectly plastic, and further, it is assumed that the incremental change in plastic true strain induced 
by an increase in true stress is proportional to the total stress deviator, rather than the incremental 
change in the stress deviator. The elastic strains are neglected in the Levy-Mises equation, and 
thus the Levy-Mises equations are only applicable when large plastic deformations occur such that 
the elastic deformations can be entirely neglected. The Prandlt-Reuss equations include both elastic 
and plastic deformations. In the present study the elastic deformations could not be neglected and 
hence the Prandtl-Reuss equations were applied, as further described below. 

Reuss assumed that the plastic true strain-increment, at any instant of loading, is 
proportional to the instantaneous stress deviation such that: 

P P P 
dei de 2 de 3 . 

— = — = — = dX (2.15) 

s'i s 2 s' 3 

P 

where dei are the plastic principal true strain increments, s\ are the deviatoric true stress 
components, and dX is a non-negative constant of proportionality which is not a material constant 
but rather varies throughout the stress/strain history. 
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Equation (2.15) implies that the incremental plastic strains depend on the total current 
deviatoric stresses and not on the incremental change in deviatoric stresses. It is also implied that 
the principal axes of stress and plastic strain increments coincide. Again, the hydrostatic stress 
component assumes no role in the plastic deformation of the material. 


For a given direction the total strain increment de l is defined as the sum of the plastic deP 
and elastic de e strain increments: 


. 1 . P , e 

dcj — dcj + dcj 


(2.16) 


The plastic strain increment deP is given by Equation (2.15), and hence [Slater, 1977]: 


. P 'A ' i 3 
de; = s . dA. = ~ 

i 2 


s'j ds 
sH 


(2.17) 


where s is the effective true stress defined in Equation (2.12) and H (also known as the “plasticity 

modulus ) is the slope of the effective true stress-plastic strain curve. The elastic strain increment 
c # 

dei is dependent on both the deviatoric and hydrostatic stress components: 

_ ^ s i (l-2v) 

de ‘ = 2G + IT' ds M (2. 1 8) 

where E is the Young’s modulus, v is the Poissons ratio, and G is the shear modulus. The shear 
modulus is related to Young’s modulus and Poissons ratio: 
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E (2.19) 

2(l+v) 

Combining Equations (2.17) and (2.18) gives the total strain increments for a given direction i: 


t 3 s i ds 
dei — 

2 sH 


ds'j (l-2v) 


2G 


- + 


dShy 


( 2 . 20 ) 


Equation (2.20) is the Prartdtl-Reuss relation applicable for initially isotropic materials to describe 
post-yield stress-strain behavior. 

A method of experimentally verifying whether it is appropriate to apply the Prandtl-Reuss 
system of equations to a given material was developed by Lode (1926). Lode introduced two 
parameters, \i and v, know n as the Lode stress and plastic strain parameters respectively: 

_ 2s 2 - si - S3 

^ si - S3 

( 2 . 21 ) 


e? - e P 
*1 e 3 


In Equation (2.21), the stresses are the total principal stresses and the strains are the plastic portion 
of total principal strains defined as : 


e p = e t - 


( 2 . 22 ) 


where e l and e e are the total and elastic principal strains respectively. 



Lode showed that if a material is elastic-perfectly plastic, then the Prandtl-Reuss equations 
imply that the two parameters will be identical, i.e., p = v. The Lode parameters are a convenient 
way to test whether a given material follows the Pandtl-Reuss equations. The Lode parameters 
have been evaluated for a number of metals and metal alloys, including iron, copper, nickel, 
aluminum, and steel alloys [Morrison and Shepard (1950); Taylor and Quinney (1951)]. Typical 
experimental results are shown in Figure 2.2. As indicated, the correlation is usually not exact, 
and often a consistent deviation from predicted behavior can be identified. This deviation from 
theory is usually attributed to anisotropy of the test material and/or to experimental errors, although 
it is also likely that the material is not completely elastic-perfectly plastic. Despite these 
inconsistencies, the Prandtl-Reuss equations are considered to be one of the best theoretical 
approaches to modeling post- yield behavior currently available. 



Figure 2.2: Lode Variables Measured by Taylor and Quinney for Copper, Aluminum, and Mild 
Steel Tubes Subjected to Combined Tension and Torsion (Taken from Hill [1950] pg 

A A\ * ® 
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Post- Yield Behavior of Anisotropic Materials 


An anisotropic form of the power law was applied to polyethylene during the present 
study. The effective true stress and effective true strain for istropic materials have been defined in 
Equations (2.12) and (2.13). These definitions no longer apply for anisotropic materials. Hill 
[1950] modified these expressions for use with anisotropic materials. The effective true stress for 
an anisotropic material in terms of the true principal stresses is given by: 


s = 



F(s?-S 3 ) z + G(sysi ) 2 + H(si-S 2) 2 
F + G + H 


(2.23) 


where: 


H + G=4; 


H + F = -^r 


G + F = 


T^ 

*3 


and si, S 2 , S 3 are the principal true stresses. Because the material is anisotropic, the normal stress 
coefficients F, G, and H are functions of the tensile yield strengths in the 1, 2, and 3 directions, 

Ti, T 2 , and T 3 respectively. 


Similarly, the effective true plastic strain for an anisotropic material in terms of the true 
principal strains is given by: 



(2.24) 
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where e^, e^, and e^ are the principal plastic normal strains and (3 = (FG + GH + HF). . Note that 

if anisotropy is negligible then Equations (2.23) and (2.24) reduce to Equations (2.12) and (2.13), 
respectively. 

For an anisotropic material, the power law formulation may be used by the substitutions of 
Equations (2.23) and (2.24) into Equation (2.14). That is, the flow law for an anisotropic material 
is given by 


s = s yield + Me’i 


where s is now given by the anisotropic form of the equivalent stress (i.e., is given by Eq. 2.23), 
and e is now given by the anisotropic form of the equivalent strain (i.e., is given by Eq. 2.24). 
Note that Eq (2.25) is identical to Eq (2.14); the variables have simply been redefined to account 
for anisotropic behavior. The anisotropic form of the power law expression was applied during 
the present study, as will be described in a following section. 



CHAPTER 3 - EXPERIMENTAL PROCEDURES 


In this study thin-walled tubes of high-density polyethylene were subjected to 
monotonically increasing axial loads and internal pressures. This biaxial loading produced both 
hoop and axial stresses within the walls of the specimen. Loading was increased proportionaly at 
user-specified rates, and the yield and post-yield response of the specimen was monitored 

throughout each test. 


Due to symmetry the axial and hoop stresses were the principal stresses in all cases. The 
axial and hoop stresses induced during a given test therefore define a radial "load path" in principal 
stress space. The load paths used during the present study are summarized in Figure 3.1. A total 
of 57 seperate tests wen: planned, involving three repeated tests along each of 19 load paths. The 
load paths are oriented at five degree increments in the first quadrant of principal stress space. 


The yield behavi or of polyethylene is highly rate-dependent. As discussed in the preceding 
chapter, in most yield- and post-yield criterion it is assumed that inelastic behavior is governed by 
the octahedral shear stress. Therefore a constant octahedral shear stress rate was maintained during 
all tests in order to minimize rate effects between individual tests and load paths. Note that 
although the octahedral shear stress rate remained constant from test-to-test, the hydrostatic stress 
rate, axial loading rate, and internal pressure rate all varied considerably from one load path to the 


next. 
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Derivation of Test Matrix 

The axial load (L) was increased at a rate of ki and the internal pressure (P) was increased 
at a rate of k2. Therefore at time t: 


L = kit 


(3.1) 

*0 

II 

K 

K> 


(3.2) 

Under this loading condition, the stresses induced in the walls of the tube are: 


a«OAL-ai = £('-^+ k 2 r j 


(3.3) 

~ lwt 

OHOOP = 02 = -£~ 


(3.4) 

ORADIAL = - \ k2t = 0 


(3.5) 

with: r = average radius of tube 

c = wall thickness 



Note that the radial stress has been equated to 
walled tubes, since the radial stress is typically 

zero. This is an appropriate assumption for thin- 
an order of magnitude less than the hoop and axial 


stresses. 
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Normally, "ai '" denotes the algebraically greatest principal stress, "02" the intermediate 
principal stress, and "03" the algebraically smallest principal stress. However, for purposes of 
consistent identification this convention was ignored in this study. Instead, "ai" was used to 
denote the axial stress, and "02" was used to denote the hoop stress, even though along some load 
paths the hoop stress was in fact the algebraically greatest principal stress. 


Substituting Eqs (3.3) and (3.4) into Eqs (2.1) and (2.2), the octahedral shear stress and 
hydrostatic stress at any time t are given by: 


^oct 


‘ /Ikii 

6c „2 r 2 


+ 6k2 2 r 2 


(3.6) 


0hl ' d 'i(to - k2r ) 


(3.7) 


The slope M of any load path shown in Figure 3.1 can be expressed as: 

= (3.8) 

The combination of load rates ki and k2 which produce a desired slope M are given by: 

ki =k27tr 2 (2M-l) (3 

Finally, given a desired octahedral shear stress rate, (Xoci/t), and slope M, the required axial load 
rate and internal pressure rate are given by: 
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( 3 . 10 ) 


(tocO 

|3c , 

(2 

1 1 J 

'r N 

3+(2M-l) 2 


( 3 . 11 ) 


An octahedral shear stress rate of 258 psi/min was selected for use since the axial loading rates and 
internal pressure rates for all loading paths falls within the range of the equipment available. A test 
matrix was defined in this fashion and is presented in Table 3.1. The numerical values given in 
Table 3.1 were obtained using Eqs (3.10) and (3.11), as well as the specimen dimensions 
presented in the following section. 

Specimen Preparation 

The material selected for testing was high-density polyethylene. Driscopipe 1000 
(PE3408) SDR26 high-density polyethylene tubing was purchased in 20 foot lengths from a 
distributor in the Seattle area. Since the polyethylene tubes were produced by extrusion, 
significant residual stresses were present in the as-received tubes. Therefore, each specimen was 
subjected to a thermal annealing process prior to testing. To the authors' knowledge, there is no 
accepted procedure for annealing high-density polyethylene. An appropriate annealing cycle was 
determined by conducting several tests at different temperatures. Starting from an annealing 
temperature of 90°C, annealing temperatures were increased in 10°C increments from one test to 
the next until the melting temperature of roughly 130°C was reached. A Blue M forced air oven 
was used. The oven was preheated so as to achieve a constant annealing temperature with minimal 
thermal gradients throughout the oven chamber. The specimens were then placed in the oven and 



Tests along Load Pat h 2 not completed. 


I 



Z 

o 

Load Path 

oc 

11.43 

5.67 

3.73 
2.75 
2.15 

1.73 
1.43 
1.19 
1.00 
0 84 
070 
0.58 
0.47 
0.36 
0.27 
0.18 
0.09 
0.00 

§ 

| 

900 

*895 

885 

867 

838 

795 

734 

657 

560 

449 

328 

202 

83 

-31 

-143 

-230 

-311 

-384 

-449 

Axial Load 
Rate 
(lh/min) 

0 

5 

10 

15 

21 

27 

34 

40 

46 

51 

54 
57 

58.2 

58.4 
57.7 

56.5 

55 
53 
51 

Int. Pressure 
Rate 

(psi/min) 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

258 

S' 

I s? 
|. % 

186 

210 

236 

260 

284 

310 

331 

348 

357 

365 

357 

348 

331 

310 

284 

260 

236 

210 

186 

1 £ 
3 - 

| ^ 


Table 3 1: Axial and Internal Pressure Loading Rates Required to Perform Constant Octahedral 

Shear Stress Rate Tests 
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annealed for one hour. The oven was then turned off and the specimens were slowly furnace 
cooled over approximately 10 hrs. 

The effect of annealing temperatures was evaluated by means of a destructive sectioning 
technique. A 1 inch wide strip was removed along the length of a 12 inch tube to evaluate axial 
residual stresses . For the hoop direction, a 1-1/2 inch long "ring" specimens was cut from the 
tube, and a 1-1/2 wide section (corresponding to a 20° arc) was cut out of the circumference of this 
ring. After a period of 72 hours (to allow for any time-dependent response to develop) the strip 
and ring specimens were examined to determine whether any deformations had occurred. 

It was found that large residual stresses existed in the tubes in the as-received condition. 
Compressive residual stresses existed in the hoop direction, while axial residual stresses were 
probably compressive at the outer diameter and tensile at the inner diameter. Figures 3.2 and 3.3 
show the results of the various annealing temperatures on hoop and axial residual stresses. Note 
that at 130°C, the specimen had been extensively deformed by the annealing process, which 
indicated that the melting temperature had been reached. For the axial direction, the amount of 
deflection from the center of the strip was measured as shown in Figure 3.4. For the hoop 
direction, the angle of deflection from the center of the tube was measured as shown in Figure 
3.5. Measured deformations were plotted functions of annealing temperature, and the results are 
summarized in Figure 3.6 and 3.7. The results show that deformation of the axial strip was 
minimized using an annealing temperature of approximately 1 15°C. However, higher annealing 
temperatures resulted in a algebraic sign change (i.e., a reversal in the curvature of the strip), 
indicating that axial residual stresses may still have existed in the tubes following an anneal at 
115°C. Furthermore, the hoop deformations were not completely eliminated even at annealing 
temperatures approaching the melting temperature of polyethylene. 
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Figure 3.2 - Hoop Sections of Annealed Tubes at Various Temperatures 
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Figure 3.3 - Axial Sections of Annealed Tubes at Various Temperatures 
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ORIGINAL TUBE 
CENTER LINE 



Figure 3.4: Measurement of deformations due to axial residual stresses. (Note: A deflection was 
defined as “positive” if the strip was concave with respect to the original centerline of 

this tube.) 




Figure 3.6: Transverse Deflections Due to Axial Residual Stresses as a Function of Annealing 
Temperature 
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Figure 3.7: Angular Deflections Due to Residual Stress as a Function of Annealing Temperature 
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The annealing treatment final,, selected for use consist of a 2 hour anneal a, U« 
followed by a 10 hr furnace cooling. it is sieved lhal Ihis annealing — — ed, 
perhaps did not completely eliminate, residual stresses present in the as-received tube 

Specimens were cut from '•parent" tubes in ,5-1/2 inch lengdis and annealed. Specimen 
dimension measurements were, taken on six different sections of the annealed tubes in or er to 

, The average radius of the annealed tubes was 1.698", and the average 

determine average values. The average ramus oi 
wall thickness was 0.157". 


Axial Load and Internal Pressure Control 

Rgure 3.8 shows the fixture used to apply axial loads and internal pressures to the tubular 
specimens. The specimens were mounted to two aluminum cylinders. The upper cylinder 
extended into the specimen in order to reduce die volume of air retired to 
specimen, which deceased the energy released in the event of catastrophic specimen ■ 
pons were machined in die lower cylinder. One was die supply port for die pressunxed 
homed air cylinder and die oiher was die pleasure transducer pen for meastuing interna, pressures. 

The tube specimens were seated on a shoulder machined in each aluminum cylinder. An O- 
ring was used to insure proper pressurixation with no leakage. After seating the tube on te 
shoulder of each cylinder end, an aluminum clamp was attached at each end to f y se t e u 
and cylinder surfaces. The Camp also helped prevent slippage of die tubes during testing y 
pressing die specimen wads into die machined step of the cylinders, thus incensing the friction 

force between the mating surfaces of the tube and cylinder. 
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Axial loads were applied using an MTS electrohydraulic fatigue testing machine consisting 
of a Model 308.01 loading frame, a hydraulic cylinder and pump, a Model 442 control unit, a 
Model 410 digital function generator, and a Transducer Inc. Model WTC-FF42-CS-10K 10 kip 
load cell. A Tektronix Model 2213 60 MHz oscilliscope was used to monitor and verify the axial 
loading. Figure 3.9 illustrates the loading frame along with the load cell and specimen. 


The MTS frame was used to apply a monotonically increasing axial load at a user-specified 
rate. The required axial loading rates have been previously listed in Table 3.1. The applied load 
rate was monitored throughout each test, and was maintained to within ±5% of the desired load 

rate in all cases. 


A control system designed to produced a linearly increasing air pressure was developed 
during this study. The system is shown schematically in Figure 3.10, and a photograph of the 
assembled system is shown in Figure 11. The pressure regulation system involved two pressure 
regulators mounted in series. First, a Victor mode SR4G regulator was mounted directly on the 
stem of a compressed air bottle, and reduced the nominal air pressure from 3000 psi to 300 psi. 
This was referred to as the "knockdown" regulator. Air flow was then directed to a Watts R1 1- 
03D regulator. This second regulator was referred to as the "control" regulator. The control 
regulator outlet pressure can be adjusted from 0 to 250 psi. An angular rotation of the valve stem 
of the control regulator results in a pressure change at the outlet The measured outlet pressure was 
found to be linearly proportional to the angular position of the valve stem. Thus, by rotating the 
valve stem at a constant angular rate the outlet pressure was increased at a constant linear rate. 


The control regulator was coupled to a Slo-Syn M093-FD14 stepper motor. This motor 
was controlled by an Anaheim Automation BLB driver, which was in turn regulated by analog 
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voltage pulses from an A/D board mounted in a Macintosh II computer. A single pulse from the 
computer caused the stepper motor shaft to turn 1.8°. A pulse train of a given frequency would 
therefore cause the motor shaft to turn at a given angular frequency. Since the shaft was coupled 
to a linear pressure regulator, a given pressure rate was produced. The entire system was 
calibrated such that a specified pulse rate resulted in a known linear pressure rate. A computer 
program was written which allows the user to produce a pulse train of a desired frequency. 
Typical plots of pressure versus time recorded along load paths 5 and 10 are shown in 
Figure 3.12. 

Data Acquistion 

A total of five analog measurement devices were used during any test. Axial strains were 
measured using an MTS model number 632.25B-20 extensometer. Hoop strains were measured 
using two diametrically opposed Applied Test Systems (ATS) LVDT displacement transducers, 
with a range of 0.5 inch. Pressures applied to the specimen was measured using a Gulton GS616 
pressure transducer. Finally, axial loads were measured using the load cell mounted to the MTS 
frame, as previously described. All five analog voltage signals were monitored using a National 
Instruments NB-MIO- 16H A/D board mounted within a Macintosh II computer. The A/D board 
was controlled using "Labview", a National Instruments graphical programming language. The 
Labview package allowed scanning rates of up to 45000 samples per second and up to 8 double- 
ended input signals. One potential difficulty with the scanning process was that all 5 channels 
were not monitored at the same time but rather sequentially. A maximum time lag of 5 msec 
occurred between scantling each channel. Since the load and strain rates imposed during all tests 
were relatively low, the 5 msec lag between measurements was negligibly small, and all 5 
measurements were treated as if obtained simultaneously. 




Figure 3.8 - Specimen Clamping Fixture 
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Figure 3. 10 - Flow Chart of Pressure Regulation System 
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Figure 3.11- Pressure Regulation System 
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The scanning process created a binary data file. 


A computer program was then used to 


convert the binary data to numerical voltages. A second program converted the voltage readings to 
stross and strain values using calibration coefficients of the data measurement devices. Axial strain 
was calculated by: (1) converting the voltage reading at . - 0 sec to a corresponding gage length. 


(2) subtracting out the voltage reading at , = 0 sec from the entire se, of voltage readings to 
determine the change in voltage, and (3) multiplying the voltage readings by the calibration 
coefficient 0.0992 inches/volt and dividing by the gage length for that test. Hoop strain was 
calculated by: (1) subtracting the voltage reading for , - 0 sec from all voltage readings for each 
LVDT, (2) multiplying tie voltages for each LVDT by their corresponding calibration coefficient, 
0 28457 inches/volt and 0.29379 inchcs/volt, to determine the displacement of each LVDT , (3) 
adding the displacement forLVDTla. timet, o tire displacement of LVDT 2 at time node, ermine 
the total change in diameter of tire tube a, time t, and (4) dividing tire change in diameter by the 


average hoop diameter of 3.396 inches. The pressure and axial load at any time , was calculated by 
multiplying the voltage reading from the pressure transducer and load cell at time t by their 
respective calibration coefficients, 59.386 psi/volt and 1000 lbs/volt. The hoop and axial stresses 
in the tube were calculated using thin-walled tube theory. 


Summary of Test Procedure 

Each specimen was mounted in die MTS frame, and a small offset preload (50 lbs) was 
applied to allow the grips to set The axial load was then rezeroed. The desired axial load rate was 
entered into the MTS controller and all of the measurement devices and pressure hoses were 

attached. 
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Next, an offset pressure of approximately 7 psi was applied. Tire Labview ptogram would 
then begin the pressure ramp and data acquisition and the MTS axial load tamp would be started 
manually. Although a small time offset existed between the start of the two ramps, the error 
induced was small. An X-Y char, recorder was used to obtain an axial load versus axial strata 


hardcopy plot, allowing for real-time monitoring of events, 
bottomed out, when the maximum pressure of the system 

failed. 


The test would end when the MTS 
was reached, or when the specimen 


As mentioned, an offset pressure of approximately 7 psi existed in the specimen at the stint 
of the test. This pressure induced small strains and stresses in bod. the hoop and axial directions. 
These strains were accounted for by plotting the hoop and axial strains versus the octahedral shear 
stress, as shown in Figure 3.13. The strain values wen. duough aero since the instrumentation 
had been zeroed prior to the start of the rest The octahedral shear stress did not go duough aero 
since there was initially pressure in the system inducing hoop and axial stresses. The two curves 
linearly flt. An offset strain at time equal aero was then added to all strain values, shifting 


were 


the entire curve so that it passed through die origin. Appendix A contains plots of the experimental 
hoop and axial stresses and strains obtained for each load path. 


Strain (in/in) 
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Octahedral Shear Stress (psi) 


Hoop Strain Curve Fit: y = -8.3268c- 4 + 8.9596c- 6 x R 2 = 0.821 
Axial Strain Curve Fit: y= 1.5582c- 4 - 7.0471e- 6 x R 2 = 0.962 


Figure 3. 13 - Axial and Hoop Strain versus Octahedral Shear Stress to Determine Offset Values of 

Strains due to Pressure Offset at time = 0 sec (Note: Data taken from Load Path 17 
Test 57.) 




CHAPTER 4 - DATA REDUCTION METHODS 


Determination of Material Properties 


Recall tat load paths 1 and 19 (see Rgure 3.1) define uniaxial stress tests, in which "pure 
axial and hoop stress, respectively, were applied. Four repeated tests were completed for these 
w „ load paths. The stress and strain values recorded during these tests wete used to detemune 
Young’s modulus E and Poisson ratio v, in bod, the hoop and axial directions. 


The experimentally determined material properties were: 


Ei= 83,493 psi 
E2= 108,540 psi 
vi2= 0.39 
V2l = 0*52 


As previously noted, the subscript ”1” has l*en used to denote the axial direction, whtle the 
su^crip, "2” is used to denote the hoop direction. Previously reported values for E and v range 
from 60,000 psi to 180,000 psi and 0.31 to 0.6. respectively [Manufacturing Chemtsts 
Association, 1957; Rodriguez, 1982; Ward, 1975). Measurements obtained dunng thts study 
therefore fall within the range of previously reported measurements. 


The measured properties show that the annealed high-density polyethylene tubes were 
anisotropic. I. can be shown drat these properties should satisfy the following inverse telanonshtp 


[Jones, 1975]: 
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V12 V21 

Ei " E 2 (4- 1 ) 

Using the values of Ei, E 2 , and vi 2 listed above, the expected values of v 2 i is 0.51. This 
compares very well with the experimental value of 0.52. For consistency, the value of v 2J in 
subsequent calculations was forced to fit the inverse relationship, i.e., v 2 i was equated to 0.51. 

Recall that the tubular specimens were produced by an extrusion process which resulted in 
significant residual stresses in the as-received tubes. It is likely that the extrusion process also 
caused the pronounced anisotropy in mechanical properties. The specimens were all subjected to a 
thermal annealing treatment prior to testing, so as to relieve the residual stresses. It had been 
hoped that the annealing process would also reduce or eliminate any anisotropic material behavior. 
This result was not achieved. The annealing process minimized residual stresses, but had little 
impact upon the anisotropy of the material. 

Prediction of Yielding - Isotropic Model 

Although the test material was clearly anisotropic, it was first modeled as an isotropic 
material. This approach may be of practical interest, since it avoids the additional mathematical 
complexities associated with anisotropic constituitive models. Average values were assigned to E 
and v, 96000 psi and 0.46 respectively, and were assumed constant throughout the material. 

Yield was based upon a 0.3% strain offset on the octahedral shear stress- octahedral shear 
strain curve. Octahedral shear stress and strain were used rather than individual axial and hoop 
stress and strain components in order to determine one value for the yield point and to bypass the 
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need to average values found on the individual stress-strain curves. This approach is in contrast to 
the method used by previous researchers [Raghava et al, 1973; Ely, 1968] wherein yield is defined 
based on axial stress and hoop stress versus effective strain curves. This latter approach results in 
two measurements of yield strength, which are subsequently averaged. A 0.3% offset was used 
because it was the most common offset in previous yield studies on polymers [Raghava and 
Caddell, 1974; Raghava et al, 1973]. 

The octahedral shear stress was calculated using Equation 2.2, repeated here for 
convenience: 


'Coct = i V(°r^ 2) 2 + ( 02 - 03) 2 + ( 03 - 0 1) 2 


(4.2) 


Octahedral shear strain is a similar function of the three principal strains: 
Yoct = |V( £ l-£2) 2 + (£2-£3) 2 + (£3-£l) 2 


(4.3) 


Since no measurement of £3 was made during testing, it was necessary to determine £3 as a 
function of £1 and £ 2 . Assuming that a plane stress state existed in the thin-walled tube, £3 is 
related to the two measured strains in the following manner: 


£3= — [El+Ell 

1-V 


(4.4) 


Octahedral shear strain can now be written: 
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y™4^f^^ 2+t22) + *' e2 V (i + :v <4 ' 5) 

It can be shown [Semeliss, 1990] that the slope of the octahedral shear stress vs octahedral shear 
strain curve, assuming isotropic material behavior, is given by: 

SH = — = G ( 4 -6) 

dy 2(v+l) 


That is, the slope of the curve is independent of load rates and is simply equal to the shear modulus 
G. This simple relationship does not hold for anisotropic material behavior, as will be shown in a 
following section. 

The procedure used to identify the yield point was as follows. Three repeated tests were 
conducted along each load. path. All measurements for a given load path were combined into a 
single data set, and a fifth-order polynomial was fit to the data. The yield point was then identified 
by calculating the point of intersection between a 0.3% strain offset with a slope equal to G and the 
fifth-order polynomial. A typical result is shown in Figure 4.1. This numerical approach 
eliminated any need for estimating the intersection visually, and defined the yield point 
conveniently and consistently from one load path to the next. Since the axial stress to hoop stress 
ratio was constant for each test, the axial and hoop stresses which correspond to the yield point 
could be calculated based upon the octahedral shear stress at yield. Appendix B contains plots of 
octahedral shear stress versus octahedral shear strain for all load paths. 


Octahedral Shear Stress (psi) 
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Figure 4.1 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 3 with 
Isotropic Assumption: Axial Load Rate = 885 lb/min. Pressure Rate = 10 psi/min. 
Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 236 psi/min. 
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Prediction of Yielding - Transversely Isotropic Model 

Hooke's law for an orthotropic (or transversely isotropic) material subjected to a plane 
stress state and referenced to the "principal" material coordinate system can be written [Jones, 
1975]: 


~ei" 


"S 11 S 12 0 - 


"or 

£2 

= 

S12 S22 0 


02 

-Y12- 


on 

O 

O 
1 


-X12- 


£3 = S13O1 + S23O2 (4.8) 


where: 


Si 



Sl2= S21 = 


V12 

El 


Sl 3 = S31 = 


V13 

El 



S23 = S32 = - 


V23 

E 2 


It was assumed that the material was transversely isotropic, and hence V13 = V23, E2 = E3, and S13 
= S12. Finally, V23 was equated to the value of the average Poisson ratio in the axial and hoop 
direction, i.e., V23 = 0 46. 


The radial strain was calculated assuming plane stress conditions: 

-E1V12 (1 +V23) £ 1 ~ (E2V12 2 + E1V23) £2 
£3 Ei - E2V12 2 


(4.9) 
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Octahedral shear strain for the transversely isotropic model can now be calculated as: 


Yoct = — — | z-Vaei 2 + Beie 2 + Ce 2 2 

3(Ei - E 2 v 12 2 ) 


(4.10) 


where: 


A = 2Ei 2 (Vi 2 2 V 2 3 2 + 2Vi 2 2 V 2 3 + Vi 2 2 + Vi 2 V 2 3 + Vi 2 + 1) - 2 EiE 2 (Vi 2 3 V 2 3 
+ Vi 2 ^ + 2v j 2 2 ) + 2E 2 2 Vi 2 4 

B = 2Ej 2 (2vi 2 v 2 3 2 + 3vi 2 v 2 3 + Vj 2 + V 2 3 - 1) + 2EiE 2 (Vj 2 3v 2 3 + Vi 2 3 
- Vi 2 2 V 2 3 + 3vi 2 2 ) - 4E 2 2 Vj 2 4 

C = 2Ei 2 (v 2 3 2 -r v 2 3 +1) +2EjE 2 (v i 2 2 v 2 3 - Vj 2 2 ) + 2E 2 2 vj 2 4 

Note that retention of the anisotropic material properties has a dramatic effect on the calculated 
value of octahedral she;ir strain. This fundamental difference between isotropic and anisotropic 
models will result in a di fference in the measured yield point, as will be pointed out below. 

It can be shown [Semeliss, 1990] that the predicted slope of the octahedral shear stress 
versus octahederal shear strain curve along each load path is given by: 


EiE 2 

/ ki 2 3r 2 k 2 2 

2 V 2 ^ 

y 2jt 2 r 2 + 2 


d Y VDki 2 + Fkik 2 +Gk 2 2 


(4.11) 


where: 
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D = 


E2 2 (vi;> + l) 2 
4 r 2 n 2 


F - — [E2 2 (-V12 2 + 1 ) + E1E2 (-V12 + V12V23 + V23 - 1 )] 


r 2 

G = 4" [E2 2 ( v :i 2 2 - 2vj2 + 1 ) + 2E1E2 (V12 - V12V23 + V23 - 1 ) + 4 Ei 2 (V23 2 
+ V23 +1)| 


Unlike the isotropic case (i.e., Eq 4 . 6 ), the slope of the stress-strain curve for a transversely 
isotropic material is a f unction of load rates and therefore of load path. 

The yield point was determined for the transversely isotropic model using the same 
procedure as described for the isotropic model. However, the slope is now a function of load 
path, and the octahednil shear strain values differ. Appendix B contains plots of octahedral shear 
stress versus octahedral shear strain for all load paths. 

Prediction of Post- Yield Behavior 

Power Law Models 

Isotropic Form: The power law expression was given by Equation ( 2 . 11 ), repeated here for 
convenience: 


s — Syj e id + Me n 


( 4 . 12 ) 
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where: 



1 

[(S1-S2) 2 + (S2-S3) 2 + (S3-S1) 2 ] 2 


1 

e =-y- [(ei-e2) 2 + (e2-e3) 2 + (e3-ei) 2 ] 2 


(4.13) 


(4.14) 


are the effective true stress and effective true strain defined as a function of the principal stresses 
and strains respectively. Assuming plane stress, then s 3 - 0, and Equation (4.13) can be 

simplified: 

„ _ ,5 (4.15) 

S = [si 2 + S2 2 + s 1 S2] z 


Also by assuming the constancy of volume, Equation (4.14) can be rewritten as: 

e =[5(ei2 + e 2 2<-eie 2 )]* (416) 

The true strains used in Equation (4. 16) are the plastic portion of the total true strain defined earlier 
and repeated here: 

eW-e' <4 ' 17) 

where e l and e e are the total and elastic true strain components respectively. The total true strain e l 
are the experimentally recorded values and the elastic true strains e c are defined using the biaxial 
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form of Hooke’s law. The principal stress and strain components in Equations (4.15) and (4.16) 
correspond to the axial and hoop directions, due to loading and specimen geometry symmetry. 

Taking the logari thms of both sides of Equation (4.12) gives: 

log(s-syield) = n loge + logM (4. 1 8) 

which is in the form of a linear equation: 

y = mx + b (4.19) 

where y * log s-syield, x = loge and b = logM. Thus, a log-log plot of s-s yield versus e can be 
used to determine the constants M and n, where M is the y-intercept and n is the slope of Equation 
(4.19). 


A polynomial curve-fit was used to achieve a smooth curve after shifting the individual 
tests to the isotropic yield points. A “new” data set consisting of axial stress strain and hoop stress 
strain values was generated using the polynomial curve-fitting equation. The effective true stresses 
and effective plastic strains were then calculated using Equations (4.15) and (4.16). A log-log plot 
of s-syield versus e was performed and a linear curve-fit was used to achieve the best correlation 
coefficient. This linear curve-fitting equation provided easy calculation of the constants M and n, 
the strength coefficient and strain hardening exponent, respectively. 

Anisotropic Form: For the anisotropic case, the power law relation is defined in the same way 
as the isotropic case given by Equation (4.12). On the other hand, the effective true stress and 
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—— — --r 

=: — 

• c fi 23) and (2.24) can be simplified as. 
isotropic, then Equations (2.23) and ( 


p s ^ + Gs 2 , + H(si-S2) 


s = 


F-i-G+H 


(4.20) 


FG + GH + HF 


(4.21) 


where: 


H + G = 


1 



H + F = 


1 



h = g 

■ nf .he tensile yield strengths in the axial and 
- r and H are now functions of the tensne y 
The coefficients, F, G, and H, arc 

hoop directions. 

(4 12) (4 20) and (4.21), the constants M and n can again be fou 
Using Equations (4.12), ( • ) 

using a linear curve-fit of a log s-s yiel d versus log * curve. 
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Prandtl-Reuss Model 


The Prandtl-Reuss equation for an isotropic material was given in Equation (2.20), 
repeated here for convenience: 

s'jds ds i (l-2v) , 

aei 2 sH " 2G + E dShyd (4.22) 

where s is the isotropic form of the effective true stress defined by Equation (4.15). Equation 
(4.22) gives the incremental change in total true strain associated with an incremental increase in 
stress. Note that by definition the plastic true strains equal zero at the yield point. Therefore the 
total true strain induced by a given stress state beyond the yield point equals the sum of the elastic 
response (including any strain offset associated with the definition of yielding) plus the sum of 
incremental increases, in strain calculated using Equation (4.22). The given stress state is 
represented by the effective true stress s , the deviatoric stress , and the hydrostatic true stress 
s hyd* 


In order to apply the Prandtl-Reuss equation in practice, one must determine the material 
constants E and v, as well as the plasticity modulus H. In the present case E and v were assigned 
average values of 96000 psi and 0.46. The plasticity modulus H was determined using data 
collected along the two uniaxial stress load paths, specifically load path 1 (pure axial stress) and 
load path 19 (pure hoop stress). Four repeated tests were performed along both of these load 
paths. The resulting eight data files were combined into a single data set and a plot of effective 
plastic true strain versus effective true stress was generated, as shown in Figure 4.2. The data was 
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Figure 4.2: Plot of Effective Plastic True Strain versus Effective True Stress, Used to Determine 
the Plasticity Modulus H. 
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then curve-fit using a fifth-order polynomial such that the derivative of this equation represents the 
inverse of the plasticity modulus H: 



de p 

d! 


(4.23) 


As mentioned above and described in a preceding section, yielding was defined in terms of 
a 3% offset in octahedral shear strain, and the axial and hoop strain components associated with 
this strain offset were added to the calculated elastic strains. 

Lode Parameters: The Lode parameters, p and v, provide a simple check to determine if the 
Prandtl-Reuss equations are valid for a given material. If the equations hold, then p and v are 
equal, as previously shown in Figure 2.2. For plane stress conditions the Lode parameters reduce 
to: 


= 2s 2 - si 
K si 

3e£ 

2e^ + e p 

where in the present case the principal stresses and strains correspond to the axial and hoop 
directions. The stresses in Equation (4.26) are the total stresses, while the plastic strains were 
derived also from experimental values using Equation (4.17). The parameters p and v were 
calculated for the the entire post-yield history, and then an average value was determined for each 
particular load path. 



CHAPTER 5 - RESULTS AND DISCUSSION 


YIELD PREDICTIONS 


Isotropic Models: The two isotropic yield criterion considered in this study were the Von Mtses 
yield criterion and the Pressure-modified Von Mises yield criterion. Experimentally deternuned 
yield points are listed in Table 5 . 1 . and are compared with predictions in Figure 5.1 Note .ha, the 
shape and size of the theoretical Pressure-modified Von Mises locus is a functton of the 
compressive to tensile yield strength ratio (see Eq 2.4). Since the compressive yteld strength o 
polyethylene was no, measured during this study, i, was neccessary to use hterature reference 
values for the compressive to tensile yield stiength ratio (OT) for high-density polyethylene. Two 
values have been reported in the literature for (OT): 1.3 [Raghava and Caddell, 19731 and 1. 
[Raghava e, al, 1973; Caddell e, al, 19741. Thus, the compressive yield strength is typtcally htgher 

than the tensile yield strength. 


The comparison presented in Figure 5.1 shows tha, tire experimental measurements were 
no, well predicted by either of the two isotropic yield criterion considered. Tire theoretical loc, 
were symmetric with respect to the 45“ line, while the experimental yield locus was skewed. Tins 
discrepancy is of course due to material anisotropy. As previously stated, the isotropic analysts 
was conducted because isotropic models avoid the additional madrematical complexities assoc, a, ed 
with anisotropic constituitive models, and hence are easier to apply m practice. The results 
represented by Figure 5.1 show tha, very significant errors are introduced by the assumption o 


isotropy. 
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Table 5. 1 - Experimental Yield Stresses for the Isotropic and Transversely Isotropic Models 
Isotropic Model Transversely Isotropic Model 


Load Path 


1 

3 

4 

5 

6 



931 

0 

835 

147 

748 

200 

871 

317 

906 

421 

981 

567 

1043 

729 

1099 

923 

1283 

1283 

1209 

1439 

1024 

1463 

928 

1600 

735 

1564 

644 

1789 

424 

1572 

323 

1794 

116 

1293 

0 

1438 


0 

199 

288 

438 

576 

749 

926 

1082 

1262 

1480 

1406 

1459 

1413 

1627 

1390 

1560 

1220 

1279 





Axial Stress (psi) 
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Figure 5.1 - Comparison of Experimental Measurements and Theoretical Yield Loci for 

the Isotropic Models 



Axial Stress (psi) 
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Figure 5.2 - Comparison of Experimental Measurements and Theoretical Yield Loci for 

the Transversely Isotropic Model 



Figure 5.3 - Standard Deviation of Pressure-modified Tsai-Hill Yield Locus from Experimental 
Yield Locus as a Function of Compressive to Tensile Yield Strength Ratio (C7T) 


! 
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It was concluded on the basis of these results that the yield behavior of high-density 
polyethylene cannot be adequately modeled using isotropic yield criterion. Further, this conclusion 
is likely to be true for the general class of semi-crystalline thermoplastic polymers. 

Anisotropic Models: The two anisotropic yield criterion considered in this study were the Tsai- 
Hill criterion and the Pressure-modified Tsai-Hill criterion. Experimentally measured yield points 
based on the transversely isotropic model are listed in Table 5.1. A comparison of the 
experimental yield locus and the predicted loci of the two anisotropic models is shown in Figure 
5.2. Once again, the compressive to tensile yield strength ratio was needed to define the theoretical 
Pressure-modified Tsai-Hill yield locus (see Eq 2.6). In this case both hoop and axial strength 
ratios are required. The same values for the strength ratios mentioned above were used (i.e., C/T 
= 1.3, 1.08), and the ratios were assumed to be the same for both hoop and axial directions. 

Comparing Figs 5.1 and 5.2, it is immediately obvious that the anisotropic models predict 
the experimental behavior much better than the isotropic models. However, by inspection alone it 
is difficult to determine whether the Pressure-modified Tsai-Hill yield criterion or the Tsai-Hill 
yield criterion best fit the experimental results. A calculation based on the standard deviation 
between measurement and prediction was conducted to determine which criterion correlated best 
with experimental results. The procedure is fully described in Appendix C. Briefly, the standard 
deviation was based on the difference between the radial distance to an experimentally measured 
yield point and the radial distance to the corresponding point on the theoretical yield locus defined 
by the Pressure-modified Tsai-Hill criterion and some given C/T ratio. In essense, the procedure 
determined the particular C/T ratio which resulted in the best fit between measured and predicted 
yield response. 
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The standard deviation as a function of the C/T ratio is shown in Figure 5.3. The results 
indicate a minimum standard deviation when the strength ratio was 1.11. As discussed above, C/T 
ratios for high density polyethylene equal to 1.08 and 1.3 have been reported. Hence, the C/T 
ratio deduced from the data collected during this study is in close agreement with previously 
reported values. 

To summarize, the yield response of annealed high-density polyethylene is best predicted 
using the Pressure-modified Tsai-Hill criterion, with a C/T ratio of 1.1 1. This criterion adequately 
models both the anisotropic nature of polyethylene, and also accounts for the effects of hydrostatic 
stress on yielding. 

POST- YIELD PREDICTIONS 

Power Law Model 

Isotropic Form: Plot s of log s-s yield versus log e curves for load paths 1 and 19 of the isotropic 
case are shown in Figures 5.4 and 5.5 respectively. The straight line is the linear curve-fit of the 
experimental data points. The equation shown in the figures is the equation of the straight line and 
the corresponding correlation coefficent R 2 . The correlation coefficient for all load paths achieved 
a value of 0.99 or bette r. 

Figures 5.6 and 5.7 show the experimentally obtained values of M and n for all load paths. 
As indicated, the values of M and n varied substantially with load path. These large variations 
imply that the plastic flow of polyethylene cannot be adequately modeled using the isotropic power 
law, since M and n must be independent of load path. The values of the strength coefficient M 
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Figure 5.6: Strength Coefficient M versus Load Path for Isotropic Power Law Model 
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Figure 5.7: Strain Hardening Exponent n versus Load Path for Isotropic Power Law Model 
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range from about 15,000 psi to over 40,000 psi , while the values of the strain hardening exponent 
n range from about 0. 36 to about 0.46. 

Since both M ;ind n varied substantially with load path, the isotropic form of the power law 
was not further pursued during the study. 

Anisotropic Form: A plot of the experimentally determined values of M and n versus load path 
for the anisotropic power law are shown in Figures 5.8 and 5.9 respectively. Once again, M and n 
do not have the same values for all load paths, although the deviations in M and n are less than in 
the corresponding isotropic analysis. Values range from M = 4500 psi to 9000 psi, and n = 0.38 
to 0.42. Although the anisotropic power law showed improvement over the isotropic version, the 
agreement between load paths was still considered poor. In particular, the deviation in M was 
considered excessive. As discussed in Chapter 2, the power law is essentially an empirical 
relationship, and even for initially isotropic materials deviation from the power law are common. 
The lack of agreement in the present study is therefore not suprising. 

Since both M and n varied substantially with load path, the anisotropic form of the power 
law was not further pursued during the study. 

Prandtl-Reuss Model 

Calculation of the Lode Parameters: An evaluation of whether the Prandtl-Reuss relations 
are appropriate for use with polyethylene was first performed by calculating the Lode parameters 
for each load path. As discussed in Chapter 4, a material is suited for use with the Prandtl-Reuss 
equations if the Lode parameters are equal. The Lode parameters calculated during this study are 
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Figure 5.11: Measured vs Predicted Plastic Strains for Load Path 1 
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Figure 5.12: Measured vs Predicted Plastic Strains for Load Path 3 
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Figure 5.14: Measured vs Predicted Plastic Strains for Load Path 5 
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Figure 5.15: Measured vs Predicted Plastic Strains for Load Path 6 
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Figure 5.17: Measured vs Predicted Plastic Strains for Load Path 8 
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Figure 5.18: Measured vs Predicted Plastic Strains for Load Path 9 
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Figure 5.19: Measured vs Predicted Plastic Strains for Load Path 10 
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Figure 5 20: Measured vs Predicted Plastic Strains for Load Path 1 1 
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Figure 5.21: Measured vs Predicted Plastic Strains for Load Path 12 
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Figure 5.23: Measured vs Predicted Plastic Strains for Load Path 14 
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Figure 5.25: Measured vs Predicted Plastic Strains for Load Path 16 
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Figure .*>.26: Measured vs Predicted Plastic Strains for Load Path 17 
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Figure 5.27: Measured vs Predicted Plastic Strains for Load Path 18 
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compared with results measured during three repeated tests (exceptions are Figs 5.11 and 5.28, 
which show results for four repeated tests along load paths 1 and 19, respectively). Finally, note 
that all strains are plotted as positive values. The actual algebraic sign in each case is noted in the 
figure caption. 

The comparison between measurement and prediction is considered suprisingly good, since 
these predictions are based on the Prandtl-Reuss equations, which are strictly valid only for 
initially isotropic materials. Nevertheless, in most cases the predicted response falls within the data 
scatter. The greatest errors are the hoop strain predictions for the first few load paths (i.e., part b 
in Figure 5.1 1 through 5.15, corresponding to load paths 1-6). In these cases the predicted hoop 
strains at high stress levels fall below the measured hoop strains. It is believed that these 
discrepancies are due to necking of the specimen. That is, for near-uniaxial loadings and at high 
stresses the specimen began to neck, resulting in an unstable increase in the measured hoop strain. 
The Prandtl-Reuss equations do not account for the necking phenomena, and hence the predicted 
hoop strains are lower than measured in these cases. 

Discussion of Pre- and Post-Yield Behavior 

On the basis of the results discussed above, it has been concluded that the pre-yield behavior of 
high-density polyethylene in extruded tubular form is anisotropic (i.e., transversely isotropic), and 
yet the post-yield behavior is effectively isotropic. Two hypotheses which may explain this 
behavior have been considered. In the first hypothesis the polymer is visualized as consisting of 
stiff crystalline regions, completely surrounded by relatively flexible amorphous material. Upon 
the application of elastic stress levels, motions/distortions occur at the molecular level in both the 
amorphous and crystalline regions of the polymer. Since the crystalline regions were oriented 
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during the extrusion process, the initial elastic response is anisotropic. Now, as the load is 
increased and plastic stress levels are induced, the molecular motions in the crystalline regions 
reach a maximum. All additional strains are associated with molecular motions which occur within 
the amorphous regions alone. Hence, the post-yield behavior is isotropic. If this hypothesis is 
valid, then no increase in the degree of crystallinity would occur as a result of the strain levels 
imposed, although the crystalline regions may become reoriented. Orientation of the crystalline 
regions due to high levels of stress (or strain) certainly occurs for polyethylene; the initial 
orientation of the as-received tubular specimens was induced during the extrusion process, for 
example. However, the strains levels induced during this study (which in all cases were < 20%) 
are relatively modest compared to those encountered in the extrusion process (in which strains on 
the older of several hundred percent are typically imposed). 

In the second hypothesis the elastic response is again due to molecular motions within 
amorphous and crystalline regions of the polymer, and hence the elastic response is anisotropic . 
However, as plastic stress levels are reached additional stress- (cm- strain-) induced crystallization of 
the polyethylene molecular structure occurs, such that initially amorphous regions become 
crystallized. In this second hypothesis the crystalline regions are oriented along the prevailing 
maximum stress or strain direction, and hence the post-yield response is identical for all load paths. 
If this second hypothesis is valid, then the degree of crystallinity would increase as a result of 
plastic deformation. Stress-induced crystallization of polyethylene is well documeted [Rodriguez, 
1989, pg 299]. However, the stress- and strain-levels associated with these documented cases are 
far higher than those imposed in this study. Thus, it is unclear whether this mechanism can 
account for the isotropic post-yield behavior observed during the present study. 


CHAPTER 6: SUMMARY AND CONCLUSIONS 


This study was devoted to an experimental and analytical evaluation of the yield and post- 
yield behavior of high-density polyethylene, a semi-crystalline thermoplastic. Polyethylene was 
selected for study because it is very inexpensive and readily available in the form of thin-walled 
tubes. Thin-walled tubular specimens were subjected to axial loads and internal pressures, such 
that the specimens were subjected to a known biaxial loading. A constant octahederal shear stress 
rate was imposed during all tests. The measured yield and post-yield behavior was compared with 
predictions based on both isotropic and anisotropic models. Of particular interest was whether 
inelastic behavior was sensitive to the hydrostatic stress level. 


The major achievements and conclusions reached during the study are itemized below. 

(1) A suitable testing facility was designed, fabricated, assembled, and fully calibarated during 
the course of the study. One of the unique features of the testing apparatus is a pressure control 
system which allows the internal pressure to be increased at a user-defined rate, such that the 
octahedral shear stress induced in the tube wall is increased at a user-defined rate. A data 
aquisition system utilizing Labview software package and Macintosh II computer was developed 
and used to collect data automatically. The system is fully operational, and will be used during the 
coming year to study the mechanical behavior of other engineering thermoplastic polymers. 

(2) An annealing process for use with extruded high-density polyethylene tubes was 
developed. Experimental measurements show that the annealing process reduces the residual 
stresses induced during the extrusion process to minimal levels. However, the annealed 
polyethylene specimens were found to exhibit a pronounced anisotropy. A 30% difference in 
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Young's modulus for the a<ial and hoop directions existed in the annealed specimens, for example. 
Thus, the annealing process removed residual stresses, but did not remove or minimize anisotropic 
material behavior. The anisotropic elastic behavior of high-density polyethylene is clearly due to 
the semi-crystalline nature of molecular structure of this polymer. 

(3) An effort was made to determine whether relatively simple isotropic yield criterion could be 
used to predict yielding to within engineering accuracies. Predictions obtained using the Von 
Mises and Pressure- Modified Von Mises criteria were compared with the measured data. Neither 
of these models adequately predicted the experimental results. It is likely that this conclusion will 
hold for other semi-crystalline thermoplastics of current structural interest (e.g., 
polyetheretherketone or polyamide), which can also exhibit a pronounced anisotropy. 

(4) The tubular specimens were subsequently treated as transversely isotropic materials, with 
the principal material coordinate system defined by the axial and hoop directions. Predictions 
obtained using two anisotropic yield criteria (the Tsai-Hill and Pressure-modified Tsai-Hill criteria) 
were then compared with the measured data. Yielding was best predicted by the Pressure- 
Modified Tsai-Hill criterion, which indicates that the yield behavior of high-density polyethylene is 
both anisotropic and sensitive to the hydrostatic stress state. It was deduced that the compressive- 
to-tensile yield strength ratio of high density polyethylene is 1.1 1, which compares very well with 
previously reported measurements. 

(5) An attempt to model the post-yield behavior using the power law were made. Both an 
isotropic form as well as an anisotrpic form was used. Neither power law formulation adequately 
modeled the post-yield behavior, in both cases the strength coefficient M and the strain hardening 
exponent n varied substantially, depending on load path. 
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(6) The post-yield behavior was best predicted using the isotropic form of the Prandtl-Reuss 
equations. This was an unexpected and suprising result, since the initiation of yielding was 
sensitive to both the anisotropic nature of the polyethylene tubes and the hydrostatic stress level. 
Neither of these factors are included in the Prandtl-Reuss formulation. The post-yield behavior of 
polyethylene is essentially isotropic, even though the initial elastic response is anisotropic. 

Two possible explanations of the isotropic post-yield response have been proposed. In the 
first hypothesis it is suggested that molecular motions associated with plastic strains occur 
primarily within the amorphous regions of the polymer, and hence the post-yield behavior is 
effectively isotropic. If this hypothesis is valid, then no increase in the degree of crystallinity 
would occur as a result of plastic deformation. 

In the second hypothesis it is suggested that stress-induced crystallization of the amorphous 
regions of the molecular structure occurs, and further that during the loading process the crystalline 
regions become oriented in the prevailing direction of maximum stress or strain. Therefore the 
post-yield response would be identical for all load paths considered during this study. If this 
second hypothesis is valid, then the degree of crystallinity would increase as a result of plastic 
deformation. Monitoring the degree of crystallinity before and after testing, as well as the 
orientation of the crystalline regions, may provide a means of determining whether either of these 
hypotheses is valid. 



LIST OF REFERENCES 


Ainbinder, S.B., Laka, M.G., Maiors, I.Y., 1965, "Effect of Hydrostatic Pressure on Mechanical 
Properties of Polymers," Polymer Mechanics, Volume 1(1), pp. 50-67. 

Bauwens, J.C., 1970, ' Yield Condition and Propogation of Luders’ Lines in Tension-Torsion," 
Journal of Polymer Science, Volume A2(8), pp. 893-901. 

Betten, J., Frosch, H.G., and Borrmann, M., 1982, "Pressure-dependent Yield Behavior of 
Metals and Polymers, " Materials Science and Engineering, Volume 56, pp. 233-246. 

Bowden, P.B., and Jukes, J.A., 1968, "The Plastic Yield Behavior of Polymethylmethacrylate," 
Journal of Materials Science, Volume 3, pp. 183-190. 

Bowden, P.B., and Jukes, J.A., 1972, "The Plastic Flow of Isotropic Polymers," Journal of 
Materials Science, Volume 7, pp. 52-63. 

Caddell, R.M., Raghava, R.S., and Atkins, A.G., 1973, "A Yield Criterion for Anisotropic and 
Pressure Dependent Solids such as Oriented Polymers, Journal of Materials Science , 
Volume 8, pp. 1641-1646. 

Caddell, R.M., Raghava, R.S., and Atkins, A.G., 1974, "Pressure Dependent Yield Criteria for 
Polymers," Materials Science and Engineering, Volume 13, pp. 1 13-120. 

Caddell, R.M., and Woodliff, A.R., 1977, "Macroscopic Yielding of Oriented Polymers," Journal 
of Materials Science, Volume 12, pp. 2028-2036. 

Carapellucci, L.M., and Yee, A.F., 1986, "The Biaxial Deformation and Yield Behavior of 
Bisphenol-A Polycarbonate: Effect of Anisotropy," Polymer Engineering and Science, 
Volume 26(13), pp. 920-930. 

Chakrabarty, J., 1987, Theory of Plasticity . McGraw-Hill Inc, p. 60. 

Christiansen, A.W., Baer, E., and Radcliffe, S.V., 1971, "The Mechanical Behavior of Polymers 
under High Pre ssure," Philosophical Magazine, Volume 24, pp. 451-467. 

Dieter, G.E., 1986, Mechanical Metallurgy . 3rd Ed, McGraw-Hill, pp 87-89. 

Ely, R.E., 1967, "Biaxial Stress Testing of Acrylic Tube Specimens," Polymer Engineering and 
Science, Volume 7(1), pp. 40-44. 

Ely, R.E., 1968, "Biaxial-stress Results for Nylon," Experimental Mechanics, Volume 8(6), pp. 
267-271. 

Freire, J.L.F., and Riley, W.F., 1980, "Yield Behavior of Photoplastic Materials," Experimental 
Mechanics, Volume 20, pp. 118-125. 

Hill, R., 1950, The Mathemat ical Theory of Plasticity. Oxford University Press, p. 20. 


Hundy, B.B., and Green, A.P., 1954, "A Determination of Plastic Stress-Strain Relations", 
Journal on Mechanics & Physics of Solids, p 16. 



96 


Jones, R.M., 1975, M echanics of Composite Materials . McGraw-Hill Book Co., pp. 37-47. 

Manufacturing Chemists Association, Inc, 1957, Technical Data on Plastics. Manufacturing 
Chemists Association Publ, p. 536. 

Mears, D.R., Pae, K.D.„ and Sauer, J.A., 1969, "Effects of Hydrostatic Pressure on the 
Mechanical Behavior of Polyethylene and Polypropylene, ” Journal of Applied Physics, 
Volume 40(1 1), pp. 4229-4237. 

Morrison, J.L.M., and Shepard, 1950, "An Experimental Investigation of Plastic Stress-Strain 
Relations," Proc. Instn. Mech. Engrs., p 163. 

Nielsen, L.E., 1974, Mechanical Properties of Polymers and Composites-Volume 2 . Marcel 
Dekker Inc, p. 265. 

Pae, K.D., and Mears, D.R., 1968, "The Effects of High Pressure on Mechanical Behavior and 
Properties of Polytetrafluoroethylene, and Polyethylene," Journal of Polymer Science, 
Volume B6, pp. 269-273. 

Pae, K.D., 1977, "The Macroscopic Yielding Behavior of Polymers in Multiaxial Stress Fields," 
Journal of Materials Science, Volume 12, pp. 1209-1214. 

Raghava, R.S., and Caddell, R.M., 1973, "A Macroscopic Yield Criterion for Crystalline 
Polymers," International Journal of Mechanical Sciences, Volume 15, pp. 967 -97 4. 

Raghava, R.S., Caddell, R.M., and Yeh, G.S.Y., 1973, "The Macroscopic Yield Behavior of 
Polymers," Journal of Materials Science, Volume 8, pp. 225-232. 

Raghava, R.S., and Caddell, R.M., 1974, "Yield Locus Studies of Oriented Polycarbonate: An 
Anisotropic and Pressure-dependent Solid," International Journal of Mechanical Sciences, 
Volume 16, pp 789-799. 

Rider, J.G., and Hargreaves, E., 1969, "Yielding of Oriented Poly(vinylchloride)," Journal of 
Polymer Science, Volume A2(7), pp. 829-844. 

Rodriguez, F., 1982, Principles of Polymer Systems . 3rd edition, Hemisphere Publishing Co. 

Sardar, D., Radcliffe, S.V., Baer, E., 1968, "Effects of High Hydrostatic Pressure on the 
Mechanical Behavior of a Crystalline Polymer - Polyoxymethylene," Polymer Engineering 
and Science, Volume 8(4), pp. 290-301. 

Semeliss, M.A., 1990, The Yield Behavior of High-Densitv Polyethylene. M.S. Thesis, 
Department of Mechanical Engineering, University of Washington. 

Shinozaki, D., and Gioves, G.W., 1973, "The Plastic Deformation of Oriented Polypropylene: 
Tensile and Compressive Yield Criteria," Journal of Materials Science, Volume 8, pp. 71- 
78. 


Slater, R.A.C., 1977. Engineering Plasticity - Theory and Application to Metal Forming 
Processes . John Wiley & Sons, pp 90-117. 



97 


Sokoinikoff, I.S., and Sokoinikoff, E.S., 1941, Higher Mathematics for Engineers an d 
Physicists . McGraw-Hill Book Co., p 136. 


Spitzig, W.A., Sober, R.J., and Richmond, O., 1976, "The Effect of Hydrostatic Pressure on the 
Deformation Behavior of Maraging and HY-80 Steels and Its Implications for Plasticity 
Theory," Metallurgical Transactions, Volume 7A, pp. 1703-1710. 


Spitzig, W.A., and Richmond, O., 1979, "Effect of Hydrostatic Pressure on the Deformation 
Behavior of Polyethylene and Polycarbonate in Tension and Compression, Polymer 
Engineering and Science, Volume 19(16), pp. 1129-1138. 

Stassi-D'Alia, F, 1967, "Flow and Fracture of Materials according to a New Limiting Condition of 
Yielding," Meccannica, Volume 2(3), pp. 178-195. 


Stassi-D'Alia, F, 1969, "Limiting Conditions of Yielding for Anisotropic Materials," Meccannica, 
Volume 4, pp. 349-364. 

Taylor, G.I., and Quinney, H., "The Plastic Distortion of Metals, Phil Trans Roy Society, p 323. 
Ward, I.M., 1975, Structures an d Properties of Oriented Polymers, John Wiley & Sons, p. 298. 


Whitfield, J.K., and Smith, C.W., 1972, "Characterization Studies of Potential Photoelastoplastic 
Material," Experimental Mechanics, Volume 12(2), pp. 67-74. 


Williams, J.G., 1980, Stress Analy sis of Polymers. 2nd Ed., John Wiley and Sons, p 77. 

Wong, R., 1990, The Post- Yield Behavior of Hieh-DensitV P olyeth y lene M. S. Thesis, 
Department of Mechanical Engineering, University of Washington, Seattle. 



98 


APPENDIX A - EXPERIMENTAL STRESSES AND STRAINS 



(a) Axial Strain versus Axial Stress (Hoop Stress = 0) 



(b) Hoop Strain versus Axial Stress (Hoop Stress = 0) 


Figure A.l - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path 1: Axial Load Rate = 900 lh/min, No Pressure, Octahedral Shear Stress 
Rate = 258 psi/min, and Hydrostatic Stress Rate = 186 psi/min. 
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Hoop Stress (psi) 



Axial Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (psi) 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.2 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 3: Axial Load Rate = 885 lb/min. Pressure Rate = 10 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 236 psi/min. 
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Hoop Streu (psi) 



Axial Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Sims (psi) 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.3 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 4: Axial Load Rate = 867 lb/min, Pressure Rate = 15 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 260 psi/min. 
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Hoop Stress (psi) 



Axial Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (psi) 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.4 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 5: Axial Load Rate = 838 lb/min. Pressure Rate = 21 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 284 psi/min. 
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Hoop Strew (pci) 

0 465 930 1395 



(a) Axial Strain versus Axial Stress or Hoop Stress 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.5 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 6: Axial Load Rate = 795 lb/min. Pressure Rate = 27 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 310 psi/min. 
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(a) Axial Strain versus Axial Stress or Hoop Stress 
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(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.6 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 7: Axial Load Rate = 734 lb/min. Pressure Rate = 34 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 331 psi/min. 
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Hoop Stress (psi) 



Axis! Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (psi) 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.7 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 8: Axial Load Rate = 657 lb/min. Pressure Rate = 40 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 348 psi/min. 
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Hoop Stress (psi) 

0 840 1680 2520 



(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (pci) 

0 840 1680 2520 



(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A.8 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 9: Axial Load Rate = 560 lb/min. Pressure Rate = 46 psi/min. Octahedral 
Sheai- Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 357 psi/min. 
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(a) Axial Strain versus Axial Stress or Hoop Stress 



(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A 9 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 

Path 10: Axial Load Rate = 449 lb/min, Pressure Rate = 5 1 psi/min, Octahedral 
Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 365 psi/min. 
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(b) Hoop Strain versus Axial Stress or Hoop Stress 
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(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A 1 1 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path 12: Axial Load Rate = 202 lb/min, Pressure Rate = 57 psi/min. 
Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate - 
348 psi/min. 
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Hoop Strew (psi) 



Axial Strew (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Strew (psi) 



Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A. 12 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path 13: Axial Load Rate = 83 lb/min. Pressure Rate = 58.2 psi/min, 
Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
331 psi/min. 
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Hoop Stress (psi) 



Axis! Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (psi) 
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(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A. 13 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path 14: Axial Load Rate = -31 lb/tnin, Pressure Rate = 58.4 psi/min, 
Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
310 psi/min. 
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Figure A. 14 


Hoop Sires* (psi) 



(a) Axial Strain versus Axial Stress or Hoop Stress 



(b) Hoop Strain versus Axial Stress or Hoop Stress 


- Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path 15- Axial Load Rate = -143 lb/min. Pressure Rate = 57.7 psi/min. 
Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
284 psi/min. 
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Axial Stress (psi) 

(a) Axial Strain versus Axial Stress or Hoop Stress 

Hoop Strew (psi) 
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Axial Stress (psi) 

(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A. 15 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
Path* 16: Axial Load Rate = -230 lb/min, Pressure Rate = 56.5 psi/min, 
Ocuthedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
260 psi/min. 
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(a) Axial Strain versus Axial Stress or Hoop Stress 
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(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A. 16 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
* Path 17: Axial Load Rate = -31 1 lb/min, Pressure Rate = 55 psi/min, 

Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
236 psi/min. 
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Hoop Stress (psi) 



(a) Axial Strain versus Axial Stress or Hoop Stress 


Hoop Stress (psi) 



(b) Hoop Strain versus Axial Stress or Hoop Stress 


Figure A. 17 

Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate 
210 psi/min. 
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(a) Axial Strain versus Hoop Stress (Axial Stress = 0) 



Hoop Stress (psi) 

(b) Hoop Strain versus Hoop Stress (Axial Stress = 0) 


Figure A. 18 - Axial Strain and Hoop Strain versus Axial Stress or Hoop Stress for Load 
b Path 19: Axial Load Rate = -449 lb/min. Pressure Rate = 51 psi/min, 

Octahedral Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 
1 86 psi/min. 
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appendix B . VIEED . AND 




Octahedral Shear Stram taLoad Path 1: Axial 
assure. Octahedral hi 
■ 186 > sitoin - 






(b) Transversely Isotropic Model 


Figure B.2 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 3: Axial 
Load Rate = 885 lb/min, Pressure Rate = 10 psi/min. Octahedral Shear Stress 
Rate = 258 psi/min, and Hydrostatic Stress Rate = 236 psi/min. 
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Figure B.3 
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y m - 43.490 + 3.7845e+4x * 7.7524*+5x A 2 + 1.0656e*7x A 3 - 7.1793e+7xM + 1.3877&f8x A 5 R A 2 = 0.996 


(a) Isotropic Model 


l 


y = 



44.135 + 4.0990e+4x - 9.1227e+5x A 2 + 1.3593*+7» A 3 - 99346e+7i A 4 + 2.1 160e<-8* A 5 


R A 2 = 0.996 


(b) Transversely Isotropic Model 


Figure B.4 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 5: Axial 
Load Rate = 838 lb/min. Pressure Rate = 21 psi/min, Octahedral Shear Stress 
Rate = 258 psi/min, and Hydrostatic Stress Rate = 284 psi/min. 




Octth'vfc-* 1 Shear Stress (psi) 


120 


y 


Figure B.5 



(a) Isotropic Model 



= - 50.039 ♦ 4.2889*4* - 1.0196*6**2 ♦ 1.6451*7**3 - 14411*8**4 + 5.0549*8**5 R*2 = 0.999 

(b) Transversely Isotropic Model 


tahedral Shear Stress versus Octahedral Shear Strain for I -oad Part & Axial 
id Rate = 795 lb/min, Pressure Rate = 27 psi/min, Octahedral Shear Stress 
e = 258 psi/min, and Hydrostatic Stress Rate - 310 psi/min. 
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Figure B.6 



. 41.225 + 3.8721e+4x - 7.8374e*5x*2 + 1.1329e+7x*3 - 9 . 4552^-7 xM + 3.2403^8x^5 R*2 = 0.997 

(a) Isotropic Model 



y = . 41.515 + 4.1926e+4x - 9.1830e+5x'2 + 1.4282e+7x / '3 - 1.2715e<-8xM + 4.6022e+8x A 5 R A 2 - 0.997 

(b) Transversely Isotropic Model 


Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 7 : Axial 
Load Rate = 734 lb/min. Pressure Rate = 34 psi/min. Octahedral Shear Stress 
Rate = 258 psi/min, and Hydrostatic Stress Rate = 331 psi/min. 
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Ocuhedral Shear Strain (0.01 in/in) 

y = - 32.901 + 3.792 5e+4x - 6.7272*+5x A 2 + 7.79l6e+«x A 3 - 4. 6830^-7 xM + 8.5759*-7x A 5 R A 2 = 0.994 

(a) Isotropic Model 



Octahedral Shear Strain (0.01 in/in) 

y = - 32.372 + 4.0716e+4x - 7.6255©f5x A 2 + 9.1073e+*x*3 - 3.3717«*7xM + 7.2984e+7x A 5 R A 2 = 0.994 

(b) Transversely Isotropic Model 


Fieure B.7 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 8: Axial 
Load Rate = 657 lb/min, Pressure Rate = 40 psi/min. Octahedral Shear Stress 
Rate = 258 psi/min, and Hydrostatic Stress Rate = 348 psi/min. 
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(b) Transversely Isotropic Model 


Figure B.8 shear Stress 

Rate = 258 psi/min, and Hydrostatic Stress Rate - 357 psi/min. 



124 



y = - 37.961 + 4.8201<*4x - 1.1612^6x^2 + 1.7207e + 7x''3 - 1.2873^8xM + 3.3491^8x^5 R A 2 = 0.993 

(a) Isotropic Model 


y = 



38.209 + 5.203 8e+4x - 1.3628**6x A 2 + Z227le+7x A 3 - 1.9105*t8xM + 6.4l05e+8x A 5 R A 2 = 0.993 


(b) Transversely Isotropic Model 


Figure B 9 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 10. 
8UIC Axial Load Rate = 449 lb/min. Pressure Rate = 5 1 psi/mm. Octahedral Shear 
Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 365 psi/min. 




(b) Transversely Isotropic Model 


Figure B.10 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 11: 

Axial Load Rate = 328 lb/min. Pressure Rate = 54 psi/min, Octahedral Shear 
Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 357 psi/min. 
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Octahedral Shear Strain <0.01 in/in) 

y = 72.185 + 3.3978e+4x - 4.5458<*5x A 2 + 1.0988e+6x A 3 + 4 0532e+7x A 4 - 3.4073»f8x A 5 R A 2 = 0.914 


(a) Isotropic Model 



0 0.5 1.0 1.5 20 25 3.0 3.5 4.0 4.5 5.0 

Ocubufral Shear Strain (0.01 m/in) 


y = 31 .096 + 4.6417e+4i - 1.4640©t6x A 2 + 3.8177»+7i A 3 - 5.7979 m-8xM + 3.4869e+9x A 5 R A 2 = 0.925 

(b) Transversely Isotropic Model 


Figure B.l 1 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 12: 

Axial Load Rate = 202 lb/min. Pressure Rate = 57 psi/min, Octahedral Shear 
Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 348 psi/min. 




(b) Transversely Isotropic Model 


Figure B. 12 Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 13: 

Axial Load Rate = 83 lb/min. Pressure Rate = 58.2 psi/min, Octahedral Shear 
Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 331 psi/min. 
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Octahwfral Shear Strain (0.01 in/in) 

. , . g.0947 + 5 0425e44* - 1.4770*46**2 4 2.6568*47**3 - 1.2050e48**4 - 1.5723*49**5 R*2 = 0.962 


(a) Isotropic Model 



y = - 7.0159 4 5.1298*44* - 1.5435e46**2 4 2.9044*47**3 - 1.5444*48**4 - 1.5222*49**5 R*2 = 0.961 

(b) Transversely Isotropic Model 


Figure B 13 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 14. 
F,gure B.13 . 3 , lb/mijl , Ra tt = 58.4 psi/nM, Octahedral 

Shear Stress Rate = 258 psi/min, and Hydrostanc Stress Rate - 310 psi/min. 
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Figure B.14 
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i'.;.176 + 7.8122e+4x - 3.9504 c-h6x*2 

(a) Isotropic Model 


+ ■ 2.1199e+9x*4 + 1.3823e+10x A 5 R*2 = 0.982 



“SSSU . ** • «■« 

. 52.788 + 7.8802e+4x - 4.0288e-M>x A 2 + 1.3156e+8x i~a> 


(b) Transversely Isotropic Model 


ihedral Shear Stress versus ^tahedr^ Shw Sttanjw L<^ ^ 
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y = 56.055 + 3.61776+4* - 3.40276+4*62 - 3.74916+7**3 + 1.10806+9**4 - 9.8166e+9**5 R*2 = 0.942 

(a) Isotropic Model 



y = 56.334 + 3.58336+4* - 3933.4**2 - 3.85556+7**3 + 1.12386+9**4 - 9.8917e+9**5 R*2 = 0.943 

(b) Transversely Isotropic Model 


•i cnirp R 15 Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 16. 
•igurc B.15 230 ib/min, Pressure Rale - 56.5 psi/mn. Octahedral 

Shear Stress Rate = 258 psi/min, and Hydrostatic Stress Rate - 260 pst/mtn. 
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y = 23.745 + 6.8111e+4* - 3.5633e+6* A 2 + 1.3294e+8* A 3 - 2.5915e+9x A 4 + 1.9514<+10x A 5 R A 2 = 0.981 

(a) Isotropic Model 



y = - 15.651 + 6.7685e*4x - 3.141fetfe*2 + 1.0368e+*x A 3 - 1.8344ef9x A 4 + 1.2705*+10x A 5 R A 2 = 0.988 

(b) Transversely Isotropic Model 


Figure B 16 - C>ctahedral Shear Stress versus Octahedral Shear Strain for Load Path 17: 

° Axial Load Rate = -311 lb/min. Pressure Rate = 55 psi/min. Octahedral Shear 

Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 236 psi/min. 





(b)Transversely Isotropic Model 


Figure B 17 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 18: 

Axial Load Rate = -384 lb/min. Pressure Rate = 53 psi/min, Octahedral Shear 
Siress Rate = 258 psi/min, and Hydrostatic Stress Rate = 210 psi/min. 




ield Stress = 678 psi 



(a) Isotropic Model 


mnn 



Ocuheckii Shear Strain (0.01 in/in) 

y ~ * 50. 8*4 + 6.47i8e+4x - Z7495e+*x A 2+ 8.0089e+7x A 3 - 1.2196e+9x*4 + 7.2970e+9x A 5 R A 2 = 0.987 

(b) Transversely Isotropic Model 


Figure B.18 - Octahedral Shear Stress versus Octahedral Shear Strain for Load Path 19: 

Axial Load Rate = -449 lb/min. Pressure Rate = 51 psi/min, Octahedral Shear 
Stress Rate = 258 psi/min, and Hydrostatic Stress Rate = 186 psi/min. 
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APPENDIX C - STANDARD DEVIATION CALCULATION 


The Pressure-modified Tsai-Hill is given by: 


H(CT1 - 0 2 ) 2 + F(C2 - O3) 2 + G(G3 - Oi) 2 + Ki<Ji + K 2 0 2 + K3O3 =1 (C. 1 ) 


with H + G = ~ 

CiM 

h + f = f V 

C2T2 

G + F = n rp 
C3T3 

K C 1 - Ti 

Kl ” C1T1 

Kl _ C2-T2 

K2_ c 2 t 2 

v, C3 ~ T3 
K3_ C3T3 


testing was completed under plane stress conditions, 03 = 0, therefore: 


(H + G) cq 2 + (H + F) G2 2 - 2Hai(J2 + Ki<Ji + K 2 o 2 = 1 (C. 2 ) 


Substituting the coefficients defined in Equation (C.l) into Equation (C. 2 ) results in: 


oi 2 cr 2 2 ai<T2 
CiTi + C 2 T 2 ' C1T1 


+ 


Cl-Ti 

C1T1 


Oi + 


C2-T 2 

C 2 T 2 


0 2 = 1 


(C. 3 ) 


It is assumed that the material exhibits the same compressive strength to tensile strength 
ratio (S) in the 1 - and 2 - directions. That is, it is assumed that: 


S = 


Ci 

Ti 


C2 

T 2 


(C. 4 ) 


Under this assumption, Equation (C. 3 ) reduces to: 
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d 2 02 2 qiQ2 S-l fc 1 Q2"\ , 

ST,2 S" 2 2 ' ST, 2 s (J, + T 2 J-‘ 


(C.5) 


For a given load j>ath of slope M, the following substitution can be made: 


ai = M 0;> 


(C.6) 


Therefore: 


M 2 ^ 2 0 2 2 Mo 2 2 S-l fM<J2 02^1 

STi 2 + ST 2 2 ‘ STj 2+ S t Ti + T 2 J _1 


Rearranging the p revious equation results in the following: 


o 2 2 f M 2 1 

s [t! 2 + T 2 2 



(C.7) 


(C.8) 


If the C/T ratio (S) and slope of the load path (M) are given, one can use the quadratic 
equation to solve for the values of hoop stress (a 2 ) which will satisfy the above 
relationship. The values of hoop stress which are the roots of the above equation 
correspond to the intersections of the Pressure-modified Tsai-Hill locus with a C/T ratio 
equal to S and die load path with the slope M. In other words, the roots of the above 
equation are the yield stresses for a given load path predicted by the Pressure-modified 
Tsai-Hill. Once the predicted hoop stress value (a 2 ) is known, the predicted axial stress is 
equal to M times the hoop stress. 
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Example: 

Given: (1) Consider Load Path 4. The slope of this load path is M = 3.73. 

(2) Consider the Pressure-modified Tsai-Hill with a strength ratio of S = C/T = 
1.08. 

(3) The tensile yield strengths in the 1- and 2- directions are Ti = 1 143 psi, and 
T 2 = 1279 psi, as given in Table 2. 

Substituting these values for M, S, Ti, and T 2 into Equation (C.8) results in the following 
expression: 

7.783e' 6 c 2 2 + 2.996447e^ a 2 - 1=0 (C 9 ) 

This results in a value for hoop stress at yield of 340 psi, and therefore the predicted axial 
stress at yield is 1268 psi. 


The experimentally measured values at yield along Load Path 3 were 288 psi (a 2 ) 
and 1074 psi (oi). The standard deviation was based on the radial distance from the origin 

to the point represented by the two stress values in principal stress space. The radius is 
defined as: 


Radius = y a i 2 + a 2 2 


(C.10) 


Therefore the experimental radius and theoretical radius as predicted by the Pressure- 
modified Tsai-Hill with a C/T ratio equal to 1.08 are: 
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Radiusexp : = 1112 
RadiuSpred = 1312 


The deviation from experimental is then: 


Deviation == 


RadiuSexp - RadiuSpred 
Radiusexp 


CC.11) 


For this example, t he deviation would be equal to -18%. The deviation of predicted from 
experimental yield stress was calculated for all 19 load paths. The standard deviation of the 
predicted yield locus based on the Pressure-modified Tsai-Hill from the experimental yield 
locus was calculated using the following relationship: 


STD = 


V 


^Deviation 2 ) 


n - 1 


(C.12) 


where n is the number of load paths. The variation of the standard deviation with the C/T 
ratio is summarized in Figure 5.3. 



